Reading the Electron Clock 
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If electron zitterbewegung is a real effect, it should generate an electric dipole field oscillating with 
the zitterbewegung frequency 2m ^c? /h. The possibility of detecting it as a resonance in electron 
channeling is analyzed. 
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INTRODUCTION 

In his seminal contribution to quantum mechanics, 
Louis de Broglie [1] conjectured that the electron has 
an internal clock oscillating with frequency vb — meC^ /h 
determined by its rest mass m^. He then deduced that 
the frequency of a moving electron observed in a lab- 
oratory will be vl = vb /t, where 7 is the relativistic 
time dilation factor. However, this idea of a particle 
clock was soon forgotten after the creation of wave me- 
chanics, wherein the de Broglie frequency ujb = Zttz^b is 
interpreted exclusively as the frequency of a wave. Be- 
sides, detecting ticks in an electron clock seemed out of 
the question because of the ultra high frequency ujb — 
7.7634 X lO^'^s^^. The possibility of direct detection as a 
resonance in electron channeling was recognized only re- 
cently and tested with positive results in an exploratory 
experiment [2]. 

When a beam of electrons is channeled along a crystal 
axis, each electron is subjected to periodic impulses from 
atoms along the axis. When the energy of the beam is 
adjusted so the period of impulses matches the period of 
the electron's clock, a resonant interaction should occur 
if the clock interacts with the atoms in some way. The 
resonant energy is easily calculated from de Broglie's hy- 
pothesis. The distance traversed during a clock period is 
d = cfi/vL = hp/{mec)'^. For the sihcon crystal used 
in the experiment, the interatomic distance along the 
< 110 > direction is d = 3.84 A, which implies a resonant 
momentum p — 80.874 MeV/c. A dip in the transmis- 
sion rate of the channeled electron beam was detected at 
81.1 MeV/c, just 0.28% from the expected resonance and 
within the estimated 0.3% limit on calibration error. 

This result may seem doubtful, as it was not antici- 
pated nor is it likely to be explained by standard quan- 
tum mechanics, despite extensive work on channeling 
theory. However, possibilities for an explanation are sug- 
gested by oscillations in solutions of the Dirac equation 
called zitterbewegung by Schroedinger [3]. Many physi- 
cists, including Dirac himself [4], have suggested that zit- 
terbewegung describes local oscillations in electron posi- 
tion that account for the electron's spin and magnetic 
moment. At first, Dirac argued for the existence of an 
electric dipole moment fluctuating with zitterbewegung 
frequency ujz = 'iujB [5], but he soon retracted that sug- 



gestion. However, the electron clock as a rotating dipole 
may be just what we need to explain the interaction pre- 
sumed in the channeling experiment, and it happens that 
ujz can account for the data at least as well as ub- In- 
deed, the two frequencies are intimately related. For a 
free particle, lub appears in the phase of a Dirac wave 
function, so observables, which are bilinear functions of 
the wave function, exhibit the frequency uJz- 

Guided by analysis of Dirac theory, a new parti- 
cle model of the electron has recently been developed 
wherein spin and zitterbewegung are complementary fea- 
tures of electron kinematics, and interaction with arbi- 
trary electromagnetic fields is included [6]. The present 
paper employs that model to explain the putative elec- 
tron clock resonance in channeling, and extend it with 
quantitative predictions for shifts in zitterbewegung fre- 
quency that can be tested with more precise experiments. 
If confirmed, the channeling experiments will thus pro- 
vide the first direct evidence that zitterbewegung is a real 
physical effect. Implications for quantum mechanics are 
likely to be profound [6] . 



THE ZITTER MODEL 

The word "zitterbewegung" is quite a mouthful, so let 
us shorten it to zitterl In developing the model of a point 
particle with zitter, geometric algebra has played an es- 
sential role, and it greatly facilitates solving the coupled 
system of equations in the model as well as comparison 
with Dirac theory [6]. However, for convenience of read- 
ers unfamiliar with geometric algebra, a synopsis of the 
model in standard tensor form is presented here (using 
natural units until calculations are in order). 

We model the electron as a point particle in spacetime 
with a lightlike history — z^{t), so its velocity 



dz^ 
dr 



satisfies 



0, 



(1) 



where r is an intrinsic time determined by the equations 
of motion (call it electron time). The particle's momen- 
tum is a timelike vector p^, and it determines a dynamical 
mass m defined by 



>0, 



and 



m 



(2) 
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The particle also has intrinsic angular momentum (or 
spin) characterized by a tensor S^^, = —Sj^^ with the 
properties 



ui'S^, = 



and 



Sjx " Sav — 0. 



(3) 



The particle is charged, so it interacts with any external 
electromagnetic field = — d^Af^. 

Particle motion is determined by a system of coupled 
equations for velocity, momentum and spin: 



du 



dr h rUf, 



q 



dr -^^^'^^ ■ 2m, 



(4) 
(5) 



Note the units with c = 1 and the two coupling constants: 
charge q and charge to mass ratio q/mg. 

The equations admit an integral of the motion. 



m = p^u^ 



me. 



2me 



(7) 



which shows that the mass is a dynamical quantity and 
specifies a shift in mass due to interaction. 



MASS, FREQUENCY AND TIME SCALING 



As suggested by this equation and confirmed by further 
analysis, the zitter vector rotates temporally with 
proper frequency cuz, so the electron's history can be de- 
scribed as a lightlikc helix winding about a timclikc curve 
determined by the momentum. The variable radius of the 
helix is inverse to the frequency. According to (7), the 
zitter frequency reduces to the 2lub for a free particle 
and corrections to this are generally small, so the radius 
of the zitter is on the order of a Compton wavelength. 

The momentum equation (5) can be put in the more 
convenient form 



dr 



= qF^^u" + 



(12) 



where the spin potential $ is given by 



2m. 
with 



qXeF"" 

n 



2me 



1 
We' 



(13) 



According to (9) the spin potential doubles as a frequency 
shift. In general the momentum will wobble about a mean 
value driven by zitter on the right side of (13), but that 
will be irrelevant to our application to experiment. 

To factor out effects of zitter on the motion, we average 
over a zitter period holding fixed. Denoting zitter 
means by an overline, we have 



Vij. = Vfj, = Up, and = 0, 



(14) 



For purposes of measurement, we must relate electron 
time to some observable time scale. A suitable relation 
is defined by the momentum vector = mw^ , which de- 
termines the velocity = dx^/dr for a timelike history 
Xfj. = x^{t), and thus associates a proper time with the 
electron motion. 

The zitter model can now be simplified and clarified 
by introducing a unit zitter vector defined by 
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where 



= — l^e + t'afiJ = -r 

n Ibjp dr 



(8) 



(9) 



can be identified as a zitter frequency and defines a zit- 
ter phase angle ip, which is closely analogous to the QM 
phase angle in the Dirac equation. 

It follows that p'^Cfj, = 0, so we can reduce the spin 
tensor to the simple form 



Moreover, the particle equation (4) is reduced to 
duu 



dr 



- = ^zeti + qF^uu". 



(10) 



(11) 



so, if Ffj.1, is not resonant with zitter fluctuations, equa- 
tion (12) reduces to 



d{mvij,) _ ^-^ ,^ 



dr 



(15) 



This is the classical equation for a structured charged 
particle subject to a Stern-Gerlach force with potential 



2me 

where Sa/s = eafftivS^v" 



(16) 



defines a spin vector with constant magnitude s^s^ = 
(7i/2)2 and s^v^" = 0. 

Finally, we can split the spin tensor into a fluctuating 
zitter component and a slowly precessing spin compo- 
nent: 



(17) 



This determines a corresponding split in the spin poten- 
tial: 



$ = qXeF^'^Vc.eff + $. 



(18) 



For physical interpretation, it is convenient to express 
$ in terms of electric and magnetic flelds E„ and 
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as seen in the rest system of the electron defined by the 
instantaneous direction of w^. Thus, 



-d • E,, 



/x.B„. (19) 

The last term on the right is the famihar Zeeman interac- 
tion energy, where fi — {q/m^)s. is the electron magnetic 
moment. The first term is the interaction energy of an 
electric dipole with d • = 0. The zitter dipole moment 
d = gAee rotates in the plane orthogonal to s with the 
zitter frequency lvz = 2m/?i. This fluctuating electric 
dipole interaction is here proposed as a physical mecha- 
nism to explain the channeling resonance. 



MOTION ALONG CRYSTAL CHANNELS 

For a quantitative account of electron channeling, in- 
teraction with the crystal is modeled by a static electric 
potential V = V{x) in the lab frame, so = — V^. 
Prom the time component of equation (12) we immedi- 
ately get the familiar energy conservation law 

E=po + V=Po + V, (20) 

where we note that the kinetic energy term must be 
equivalent to its zitter average. 

The spatial component of equation (12) can put in the 
form 

dp 



dt 



= _V(F + 7-i$), 



(21) 



where the time dilation factor 7 sd mq has been factored 
out. In the channeling domain the maximum crystal po- 
tential is a few hundred electron volts at most, so in the 
80 MeV region of interest to us, the effective electron 
mass M = jme = E/c^ is constant to an accuracy of 
10-^ and 7 = 158. 

In axial channeling electrons are trapped in orbits spi- 
raling around a crystal axis. To a first approximation, 
the crystal potential can be modeled as the potential for 
a chain of atoms, so it has the form 

F(x) = V{r, z) = U{r)P{2nz/d), (22) 

where x(t) = r + zcr^ is the particle position from the 
first atom in the chain, with r = |r|. The longitudinal 
potential P{2wz/d) = P{coot) is periodic with a tunable 
frequency loq = 2'Kz/d that varies with the energy E. 

Our problem is to calculate perturbations on the trans- 
verse component of the momentum vector, as that can 
remove electrons from stable orbits in the beam. The 
transverse component of equation (21) has the familiar 
form of a nonrelativistic equation: 

M^ = -rV' -j-\V<f )±, (23) 

where V' = drV = U' P. Ignoring the zitter perturbation 
for the time being, we seek to ascertain the effect of the 
periodic factor P{ujot) on the orbital motion. Oscillations 
in the transverse velocity can be ignored. 



Rosette channeling with atomic perturbations 

For analytic simplicity, we approximate the potential 
by the first two terms in a Fourier expansion with respect 
to the reciprocal lattice vector: 



V = U{r){l+cosLOot), 



with 



Uir) = -fcln[l + (Ca/r)^ 



(24) 



(25) 



where the first term is Lindhard's potential for a uni- 
formly charged string [7, 8], and the coefficient of the 
second term is set to make the potential vanish between 
atoms. The parameter a = 0.190A is the Fermi-Thomas 
screening radius, and the constant = 3 is a fairly ac- 
curate fit over the range of interest. For silicon {Z = 14) 
the coupling constant has the value k = Ze'^/d = 52.5 
eV. In most channeling applications the second term is 
ignored; however, its periodicity is essential for resonance 
in the zitter perturbation investigated here. 

It will be convenient to represent the radius vector in 
the complex form r = re'^ , where the imaginary i is gen- 
erator of rotations in the transverse plane. Then equation 
(23) assumes the complex form. 

= -U'Pe''^ = - — (1 + cosa;ot)r. (26) 
dt-' r 

We are interested only in radial oscillations, so we use 
conservation of angular momentum L = Mr^O to sepa- 
rate out the rotational motion. With the periodic driving 
factor omitted for the moment, equation (26) admits the 

energy integral 



El. = \Mf^ + W{r). 



where 



W{r) 



+ C/(r). 



(27) 



(28) 



2Mr2 

Let us expand this around a circular orbit of radius ro, 

and for quantitative estimates take ro = 0.50A as a 
representative intermediate radius. For ro = re*^° and 
C/q — U'{ro), equation (26) gives us 



In terms of a; = 
gives us 

1 



31.7 







(29) 



J Mro 80.9 X 10^ 
r — ro, second order expansion of (27) 



Ei_ = ^M[x^ + nlx''] + Wo, 



(30) 



where, for mass M at the expected resonance, we have 



wi 

M 



Mro 
3 X 31.7 



76.0 



80.9 X 106 



(31) 
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Thus, we have flo = 4.21 x 10^"'' s^^, which is close to 
^0 = 4.75 X 10^^ from (29) and much less than the ex- 
pected resonant frequency Wq = cje/i = 4.91 x 10^* s~^. 
Note that the distance traveled in one orbital revolution 
is dr = 27rc/0o = 3.97 x lO^A = 0.397/;im. Thus, the 
orbit makes 2.52 revolutions in passing through the one 
micron crystal, so the orbital revolutions arc only weakly 
coupled to the high frequency radial oscillations. 

Now differentiating (30) and reinserting the periodic 
driving factor, we obtain the desired equation for radial 
oscillations: 

x + nl{l + cosujot)x ^0. (32) 
This equation has solutions of the general form 



a;(t) = e'^* J2 ««« 



,inu}ot 



(33) 



Since CIq << luq the solution is well approximated by first 
order terms, which yield the particular solution 

x{t) = a{cosflt) cosLJot = ^(cosa;+t+cosa;_t),(34) 



with fl 



2 



and = LUo±n. 



(35) 



This describes a harmonic oscillator with high frequency 
Wo and a slowly varying amplitude with frequency f2, 
which is equivalent to a sum of two oscillators with fre- 
quencies Lu± separated by — lu^ = 2Cl. 

We shall see that, at the resonant frequency coq = 
sjObM, the frequency shift Q. = (0.857 x 10~^)a;o is the 
right order of magnitude to contribute to experimental 
effects. Moreover, this quantity has been estimated at 
the particular radius r = 0.50A, and it may be larger by 
an order of magnitude for smaller radii of experimental 
relevance. Accordingly, a distribution of values will 
contribute to the experiment. 



Zitter resonance 

Now we are prepared to consider the effect of zitter per- 
turbations on the orbit. Inserting the crystal electric field 

with its relativistic transformation into (19) and keeping 
only the largest term, it can be shown that [6] 

$ = -Ae7C/'Pe-f = -Ae7;7'Pcos(wzf/7 + (5).(36) 

By the way, we don't need to use the spin equation of 
motion (6) in our calculations; we only need the fact 
that it implies the unit zitter vector e rotates rapidly 
in a plane that preccsses slowly with to the spin vector 
s. We have dropped the spin contribution to because 
a more elaborate experiment would be needed to detect 
it. Also, according to equation (9), the spin potential 
(36) contributes to an oscillating shift in zitter frequency, 



but that effect is too small to be detected in the present 
experiment [6] , so in the rest of our analysis we suppose 

LOz = Re- 
inserting (36) into (23) with a convenient choice of 
phase and writing w = ojz/j, we get the equation of 
motion 

(Pr U' A 

= -r— (l + coswoi)(l + -^coswt), (37) 

where R = -U'/U" = r[l + iCa/r)^]/[3 + (Ca/r)^] is 
an effective radius. Ignoring the amplitude modulation 
as determined in (34), we can reduce this to a radial 
equation 



X + Wq(1 + cosa;t)x = 0. 
K 



(38) 



Of course, we can replace loq in this equation by lu± to 
get two separate resonant peaks. 

For h = \e/R constant, (38) is Mathieu's equation, 
for which there are well known methods of solution. So 
let us solve it using h = Xe/Ro = 1.931 x 10-^/0.208 = 
9.283 x IQ-"^ as an approximation. For small values of h 
such as this, it can be shown that (38) has resonances at 

2wo = nw, for n = 1, 2, We demonstrate that for the 

first order resonance by assuming a solution of the form 
[9] 



Lot 



cut. 



X = e [a cos — + 6 sin — 



(39) 



where the so-called Floquet exponent s must be real and 
positive. 

Looking for solutions with to = 2ujq + e close to the 
supposed first order resonance, we insert (39) into the 
differential equation (38) and keep only first order terms 
to derive the conditions 

sb-^[e-^hLOo]a = 0, sa+^[e+^hwo\b = 0.{40) 

These equations can be solved for the coefficients pro- 
vided 



(41) 



Thus, we do indeed have a resonance with width 

Aw = 2e = hujo = 9.282 x IQ-^wq, (42) 



or 



Ap = h{80.9 MeV/c) = 0.751 MeV/c. (43) 

And for the amplification factor at resonance we have 
1 



st = -hujQt = —hn = 1.46 x 10 n, 
4 2 



(44) 



where n is the number of atoms traversed in a resonant 

state. Since ln2 = 0.693, this implies that the amplitude 
is doubled in traversing about 50 atoms. 
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Now note that the value of ft, = Xe/Ro used in (42) 
and (44) applies to only the subclass of orbits for which 
vq = 0.50A. For smaller radii the values can be much 
larger. In principle, to get the width of the ensemble of 
orbits we should replace h in (42) by its average h over 
the ensemble. However, the result will probably not differ 
much from the typical value we have chosen. 

Similarly, the amplitude factor in (44) will have a 
distribution of values, and the doubling factor will be 
reached much faster for orbits with smaller radii. Pre- 
sumably, random perturbations (such as thermal fluctua- 
tions of the nuclei) will limit the resonant state coherence 
length to some mean value n. Consequently, states with 
smaller ro will be preferentially ejected from the beam. 

Even more to the point, the perturbation parameter 
h = Xe/R is not constant as we tentatively supposed but 
increases rapidly as the electron approaches a nucleus. 
In resonance the value of h close to each nucleus domi- 
nates the perturbation, so its effective mean value is much 
smaller than the estimate for constant ro . Evidently, res- 
onant interaction may eject electrons with small ro in just 
a few atomic steps. 

DISCUSSION AND CONCLUSIONS 

The predicted resonance width in (43) is in fair agree- 
ment with the width in the channeling experiment data 
[2], considering uncertainties in the value of h and such 
factors as thermal vibrations that may contribute to 
damping. Damping can only narrow the width and de- 
stroy the resonance if it is too severe. 

We need to explain how the orbital resonance is man- 
ifested in the experimental measurements. Two scintil- 
lators, SC2 and SC3, were employed to detect the trans- 
mitted electrons. The larger detector SC3, with a radius 
about 3 times that of SC2, served as a monitor while the 
smaller detector served as a counter for a central por- 
tion of the beam. The measured quantity was the ratio 
of SC2 to SC3 counts as momentum was varied in small 
steps (0.083%) over a 2% range centered at the expected 
resonance momentum 80.9 MeV. An 8% dip was observed 
at 81.1 MeV. Orbital resonances may contribute to this 
effect in at least two ways: first, and perhaps most im- 
portant, by increasing the probability of close encounter 
with a nucleus that will scatter the electron out of the 
beam; second, by increasing the duration of eccentric or- 
bits outside the central region. 

For the most part, resonant interactions will be 
strongest on electrons confined to the central region. 
However, it should be noted that in silicon each < 110 > 
channel has a nearest neighbor chain only 1.36 A away 
with the location of atoms shifted by precisely half an 
atomic step 1.92A. Considering the slow precession of 
a channeled orbit at the resonant frequency, this chain 
will resonate with it for hundreds of atomic steps. In 



fact, the interaction might lock onto the orbit to prevent 
precession during resonance. This could indeed have a 
substantial effect on resonant channeled electrons, so it 
deserves further study. These observations suffice for a 
qualitative explanation of the observed dip. A quantita- 
tive calculation will not be attempted here. Indeed, for 
a fully quantitative treatment the string potential (25) 
must be replaced by a more realistic crystal potential, 
for, as Lindhard [8] has shown, it breaks down at atomic 
distances near the screening radius a = 0.190A. Even so, 
the string approximation is useful for semiquantitative 
analysis of zitter resonance, as we have seen. 

The most problematic feature of the experiment is the 
0.226 MeV/c difference between observed and predicted 
resonance energies. If estimation of the experimental er- 
ror was overly pessimistic, that indicates a physical fre- 
quency shift. The most likely origin for such a shift is 
the frequency split in (34). The experiment was not suf- 
ficiently accurate to resolve separate peaks for the two 
frequencies, so the peaks would merge to broaden the 
measured resonance width. However, the peak for u)+ is 
likely to be higher than the peak for a;_ owing to greater 
probability for ejection from the beam. Hence, the center 
of the merged distribution will be displaced to a higher 
frequency. If this explanation is correct, then an increase 
in experimental resolution will separate the two peaks, 
and their relative heights will measure the relative prob- 
ability of ejection at the two frequencies. 

If the idea of zitter resonance is taken seriously, there 
are many opportunities for new theoretical and experi- 
mental investigations. Increasing the resolution by three 
orders of magnitude will open the door to refined studies 
of frequency shifts, line splitting, spin effects and Zee- 
man splitting, all of which are inherent in zitter theory 
[6] . As has been noted, the most straightforward predic- 
tion of the zitter model is a second order resonance near 
161.7MeV/c. In a first approximation, it can be analyzed 
in much the same way as here, though removal of elec- 
trons from the center of the beam may be enhanced by 
such processes as pair creation. 

Classical particle models have long been used for chan- 
neling calculations with considerable success. Besides be- 
ing simpler and more transparent than quantum mechan- 
ical models, they even give better results at high energies 
[7]. For present purposes, the zitter model differs from 
the usual classical model only by the zitter perturbation 
in equation (21). It is unlikely that the zitter perturba- 
tion effect can be explained by standard quantum theory, 
though it might be by a subtle modification of the Dirac 
equation [6]. The zitter model has many other physi- 
cal implications [6] , though channeling is so far the most 
direct way to test it. 

All things considered, we conclude that the quasiclas- 
sical zitter model provides a plausible explanation for the 
reported experimental data. This is ample reason to re- 
fine and repeat the experiment to confirm existence of 
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the resonance. 

I am deeply indebted to Michel Gouanere for explain- 
ing details about his channeling experiment and data 
analysis. 
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